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A model kinetic equation is solved exactly for a special stationary state describing nonlinear 
Couette flow in a low density system of inelastic spheres. The hydrodynamic fields, heat and 
momentum fluxes, and the phase space distribution function are determined explicitly. The results 
apply for conditions such that viscous heating dominates collisional cooling, including large gradients 
^ ' far from the reference homogeneous cooling state. Explicit expressions for the generalized transport 

coefficients (e.g., viscosity and thermal conductivity) are obtained as nonlinear functions of the 
coefficient of normal restitution and the shear rate. These exact results for the model kinetic 
equation are also shown to be good approximations to the corresponding state for the Boltzmann 
, equation via comparison with direct Monte Carlo simulation for the latter. 
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I. INTRODUCTION 

+-> , 

An idealized representation of rapid flow granular media is given by the Boltzmann equation for a gas of smooth, 
hard, inelastic spheres.EJ This kinetic equation provides a basis for the study of a wide range of transport properties, and 
. significant progress has been made in recent years using numerical methods such as direct Monte Carlo simulation^ In 
contrast, analytic studies have been limited to homogeneous states or to those with small spatial gradients. uu Even 
in the homogenous case there is no exact solution corresponding to the Maxwellian distribution for elastic collisions, 
1 . ■ Recently, the method of kinetic models has been proposed for practical access to transport in more complex states.cl 
This method has proven successful for elastic collisions, where several exact solutions to the model kinetic equations 
have been obtained for spatially inhomogeneous states far from equilibrium. Ut3 Such methods are potentially more 
important for granular gases since the states of interest are typically driven by external boundary conditions, posing 
l— ~~ ', intractable difficulties for solution to the Boltzmann equation. One illustration is the exact solution of a kinetic model 
' for uniform shear flow, where the rheological properties were calculated and shown to be in good agreement with tho. 

obtained from simulation of the Boltzmann equation at low densitytj and the Enskog equation at finitp-d|6iisities.l 
i Uniform shear flow is perhaps the most extensively studied inhomogeiieaus state for inelastic particles £j~I13 On the 
£""»» | other hand, a more realistic shearing state is the planar Couette AowJUlj where temperature and density gradients 
C*~) i coexist with the velocity field. The objective here is to provide an exact solution to a model kinetic equation 
corresponding to Couette flow with arbitrarily large temperature and flow field gradients, and arbitrary inelasticity. 
The results extencLa previous exact analysis of momentum and heat transport far from equilibrium for a gas with 
clastic collisions .ErtiJ Comparison with simulation of the Boltzmann equation for inelastic collisions again shows good 
agreement, confirming that the kinetic model is not only instructive but practical as well. 

As noted above, exact solutions in kinetic theory for spatially inhomogeneous states are exceedingly rare. When, 
furthermore, such a solution corresponds to a hydrodynamic state far from equilibrium a unique and important 
I ■ benchmark is obtained for both conceptual and computational issues. In the case of granular gases there are two 
interesting examples: uniform shear flow and Couette flow. Both allow controlled discussion of nonlinear rheological 
properties that are important for a wide class of real granular flows. Uniform shear flow is a useful idealization, but 
the combined heat and momentum transport of Couette flow considered here is more realistic and a stronger test of 
hydrodynamic transport. For uniform shear, the condition of stationarity imposes a relationship between the shear 
rate and the coefficient of restitution, at fixed temperature. However, for Couette flow the shear rate and the coefficient 
of restitution are independent variables since the temperature profile is allowed to change. The parameter space for 
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testing hydrodynamics is now two dimensional, with uniform shear flow recovered as a limiting line corresponding to 
zero curvature for the temperature profile. This is discussed further below. 

In light of some earlier speculation that a hydrodynamic description for granular flows might be limited to weak 
dissipation and/or weakly inhomogeneous states, the results of this paper provide an example to the contrary: a) 
a hydrodynamic description applies, since all space dependence of the heat and momentum fluxes occurs via explicit 
functionals of the hydrodynamic fields and the macroscopic balance equations become a closed set of hydrodynamic 
equations which determine the non-trivial space dependence of the hydrodynamic fields; b) the hydrodynamic descrip- 
tion applies even at strong dissipation and strong inhomogeneity (i.e., beyond the Navier-Stokes limit). The exact 
heat and momentum fluxes are characterized by generalizations of Fourier's law and Newton's viscosity where the 
thermal conductivity and viscosity are functions of the shear rate. The viscometric functions (normal stresses) are 
spatially constant, but non-trivial functions of the shear rate as well. 

These strong results are obtained in the context of a simplified model of the Boltzmann equation. The mathematical 
and physical basis for this model as a good representation of the Boltzmann equation is discussed in Ref. |6| However, 
application of the model far from equilibrium as done here raises the question of its limitations, both quantitative 
and qualitative. The Monte Carlo simulations of the Boltzmann equation provide a more secure basis for the study 
of this macroscopic state. The good agreement obtainetLin this paper for this complex hydrodynamic state extends 
that previously demonstrated for uniform shear flowjij'tij providing additional support for future applications of the 
kinetic model to address other realistic boundary driven problems (e.g., vibrated columns or vertical flow through a 
chute in a gravitational field) . Additional concerns are associated with the stability of the stationary state obtained 
here, particularly for large spatial gradients and strong dissipation. No analysis of this problem is provided here (no 
indication of instability is seen in the Monte Carlo simulation, although the imposed planar symmetry may suppress 
some possible mechanisms). 

This paper is organized as follows. In the next Section, the state of macroscopic Couette flow is first considered 
at the level of the balance equations for mass, energy, and momentum. Based on previous analysis of this problem for 
elastic collisions, it is postulated that the balance equations support a solution with constant pressure and constant 
(dimensionless) shear rate. For consistency, it is shown that the heat and momentum fluxes must be given by 
generalizations of Fourier's heat law and Newton's viscosity law, as noted above. The temperature and velocity 
profiles are then determined in terms of the coefficients in these fluxes. To confirm that this self-consistent ansatz for 
the fluxes is correct, and to determine the explicit forms for the transport coefficients, an exact solution to the kinetic 



equation for this macroscopic state is constructed in Section [II. Only the main results are quoted in this Section, 
while the detailed analysis is provided in the Appendices. The primary results of this analysis are summarized as 
follows: The heat flux is given by a generalized version of Fourier's law which is linear in the temperature gradient 
without any restriction on this gradient being small. The proportionality constant is a nonlinear function of the 
shear rate, including anisotropy effects inducing a heat flux in the direction of flow as well as in the direction of the 
temperature gradient. The momentum flux components are characterized by three scalar functions of the shear rate, 
a shear viscosity and two viscometric functions. All of these properties are applicable even far from the reference 
homogeneous state and for all values of the shear rate. The only restriction is to conditions such that viscous heating 
dominates collisional cooling so that the temperature in the bulk is higher than that at the walls. These exact results 
for the kinetic model are compared with direct Monte Carlo simulation results for the Boltzmann equation in Section 
[V. The agreement is found to be good. Finally, the results are summarized and discussed in Section M. 



II. HYDRODYNAMIC DESCRIPTION 

Consider a low density granular gas of smooth hard spheres [d = 3) or disks (d — 2) of diameter a and mass 
m. Collisions between particles are characterized through a constant coefficient of normal restitution a with values 
< a < 1 , the largest value corresponding to the elastic limit. The macroscopic balance equations for mass, energy, 
and momentum are 

D t n + riV -u = 0, (1) 



AT+^(P ij -V J -M i + V-q) =-CT, (2) 



D tUi + -nun 'V,/',, = 0, (3) 

where n is the density, T is the granular temperature, u is the flow velocity, and D t = dt + u • V is the material 
derivative. In addition, ( is the cooling rate (related to the collisional energy dissipation), P is the pressure tensor 
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(related to the transport of momentum), and q is the heat flux (related to the transport of energy). Whenever 
these fluxes can be expressed in terms of the hydrodynamic fields, Eqs. (|l|)-(§) become a closed set of hydrodynamic 
equations for these fields. 
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FIG. 1. Sketch of the geometry of the system. The velocity and temperature profiles correspond to the solution of the 
kinetic model for a coefficient of restitution a = 0.8 and a reduced shear rate a * = 0.95 [cf. Eqs. @ and @]. The units on 
the scale are such that m = 1, To = 1, and T)a(To)/p = 1, where To is the maximum temperature across the system. 



The specific problem considered in this paper is steady Couette flow [cf. Fig. [jj. The gas of inelastic hard spheres 
is enclosed between two parallel plates at y = ±L in relative motion along the s-axis and maintained, in general, at 
different temperatures. The resulting flow velocity is along the x axis and, from symmetry, it is expected that the 
hydrodynamic fields vary only in the y direction. The pressure tensor and the heat flux characterize a macroscopic 
state with combined heat and momentum transport. The objective here is to determine the hydrodynamic fields. 
Under the above conditions, Eq. (pi) is identically satisfied, while the balance equations (0) and (p|) yield 
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du x 



dq, 
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(4) 
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Equations (Q)-(|(j) are still exact. The temperature T is defined in the usual way in terms of the average kinetic 
energy, and is related to the pressure p through the ideal gas law 



p = nT = ~P U 
a 



(7) 



The second equality follows from the definition of the pressure tensor (see next Section). There are three inde- 
pendent hydrodynamic fields, which are taken here to be the pressure, temperature, and the x component of the 
flow velocity (the other d — 1 components vanish). The boundary conditions impose a global shear rate given by 
[u x (L) — u x (-L)] /2L. It is useful to introduce a dimensionless local shear rate a 



1 du r 



(8) 



Vq dy 

Here v$ is a convenient collision frequency v$ = p/rjo^ where rjo is the Navier-Stokes shear viscosity in the elastic 
limit. The variable a will be one of the control parameters for the Couette flow state. It has the simple physical 
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interpretation as the ratio of the mean free path i = vq/vq (where vq = v/2T/m is the thermal velocity) to the relevant 
hydrodynamic length h = vq(8u x / 'dy)^ 1 . We now look for special solutions to the macroscopic balance equations 
characterized by constant pressure and constant a, 

p = constant, a — constant. (9) 

To provide a closed set of equations for the hydrodynamic fields, the heat and momentum fluxes must be expressed 
in terms of these fields. Forms for these fluxes are postulated in thiSrSecjtion, as a generalization of those obtained 
from a similar analysis of Couette flow for a fluid with elastic collisions .□til Their verification is given in the following 
Section and in Appendix Momentum transport is typically characterized by three scalar rheological functions, 
the shear viscosity and two viscometric functions. The xy element for shear stresses is represented in the form of a 
generalized Newton's viscosity law 



dy 



P xy = -F v (a 1 a)r ] o^. (10) 



The function F v {a,a)rj{) is the generalized shear viscosity, where the dissipation and all nonlinear rheological effects 
are included through the factor F v (a,a). The normal stress components are specified in terms of the dimensionless 
viscometric functions 

^ 1 {a,a) = Pyv ~^ xx , * 2 (a,a) = Pzz - 2 Pm . (11) 
pa z pa z 

The heat flux is represented as 

q t = —XijVjT. (12) 

The tensor Ay- depends on two independent scalar "transport coefficients" 

Ay = A [F A (a, a)S lj + $ (a, a) a„-] , (13) 

where a%j = aSi X Sj y . The first term of Eq. ( |l3| ) yields a generalization of Fourier's law with the thermal conductivity 
\oF\(a, a) modified by the nonlinear rheological factor F\(a,a). The second term of ( |l3| ) provides information on 
the anisotropy created by the Couette flow since it gives a heat flux along the x-axis due to a thermal gradient 
parallel to the y-axis. This effect is also nonlinear with no analogue at Navier-Stokes order, being first order in 
both the shear rate and the temperature [gradient. In fact, Eqs. ( |l2|) and ( [l3|) are a generalization of the heat 
flux constitutive equation to Burnett ordeinJ for the geometry of the problem in the elastic case, according to which 
$(a = 0,a = 1) = -{8 5 - 9 i d\nrjo/d\nT)r 1 o/2mX = -3.5 (Boltzmann equation) or - 2.8 (BGK model). 

Since p and a are constant, Eqs. ( |l0| ) and (|lT|) are consistent with Eqs. (||) and (^). Dimensional analysis shows 
that 

770 cxT 1 / 2 , AqcxT 1 / 2 , (^T^ 2 , (14) 
where it is understood that the independent variables are p and T. Then Eq. (Q) becomes 

T I —T I — = ^Pr 7 , (15) 

oy oy % 



where 7 (a, a) is the dimensionless constant defined by 



Va („ 2 jp d 



2m\ F\ \ 2p ) 

and Pr = (d + 2)7y /2mAo is the Prandtl number. Equation ( p"5| ) is a closed equation that determines T = T(a, a, y). 
With this known, the velocity field is determined from the condition a = constant 

du x p 

Before giving the solution, it is instructive to express the temperature in terms of the flow velocity: T(y) — > T(u x ). 
Then Eqs. (O) and dTJ) give directly 
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d 2 T 



2m 
j-Pr7. 



(18) 



The sign of the constant 7 is a result of the competition between viscous heating and inelastic cooling (represented by 
the first and second terms, respectively, on the right-hand side of ([if])). If the dissipation is sufficiently low and/or the 
shear rate is sufficiently large so that £ < (2pF Tj /dr]o)a 2 , then 7 > and the temperature profile is convex. Equation 
( |l8| ) implies that the temperature is simply a quadratic function of the flow velocity, 



T = T a [l-c 2 (u x -u ) 2 ] 



2 _ 



(19) 



where u$ is an arbitrary constant and To is the temperature at the point where dT / du x = 0. Here and below 
attention is restricted to the case 7 > which implies that the shear rate must be larger than a certain critical value, 
a > a c (a). The relationship of T and u x can be viewed as a nonequilibrium "equation of state." The explicit 
y-dependence of u x and T can be easily obtained from Eqs. (|l9|) and (|l7|), where the constants are determined by the 
boundary conditions. For instance, suppose that y = yo is the point where the temperature presents an extrcmum 
and that, without loss of generality, u x (yo) = uq — 0. Each particular situation is then characterized by the uniform 
pressure p, the local temperature T(yo) = Tq and the shear rate a. Equation (|l7f ) becomes 



1 



2 2 N 
C U_ 



1/2 du x 



p 



dy ?7o(T)) : 



whose (implicit) solution is 



y = yo + 



2cpa 



sin 



_1 (cux) 



1 



2 2 
c u~. 



1/2 



(20) 



(21) 



Note that Eqs. ( |l9| ) and ( f2l|) are only meaningful for u 2 < c~ 2 , which implies \y — yo\ < 7r?7o(To) / 4cpa. The velocity 
profile ( ^l| ) and the temperature profile ( |l9[ ) are shown in Fig. |l] for a representative case (a = 0.8, a = 0.95) with 
yo = 0. The actual shear rate du x /dy is practically constant across most of the system, but it rapidly increases as 
the temperature becomes much smaller than To [cf. Eq. (p"7|)]. 

To summarize, the conditions of constant pressure and shear rate, Eq. (^), together with the constitutive relations 
for the heat and momentum fluxes, Eqs. (|l0|)-([i"3|) allow an exact stationary solution to the macroscopic balance 
equations. In the next Section the assumed conditions and constitutive relations are confirmed by an exact solution 
to the kinetic equation. In addition, this solution provides explicit expressions for the transport coefficients F rp ^i^j 
F\, and $ as functions of a and a. 



III. KINETIC EQUATION AND EXACT SOLUTION 



Consider a low density granular gas of smooth hard spheres in d dimensions. At sufficiently low density the phase 
space distribution function, /(r, v, t), is determined from the Boltzmann kinetic equation modified appropriately for 
inelastic collisions. Due to the mathematical complepdty of this equation, analysis has been limited to perturbative 
approximations for small spatial inhomogeneities.pTJ To describe more general nonequilibrium states it is useful 
to replace the Boltzmann equation with a more tractable model kinetic equation constructed to preserve its most 
important qualitative features (e.g., macroscopic balance equations). The model kinetic equation chosen for analysis 
here is 



dtf- 



V/ = -!/(/ - ft 



\cd v • [(v 



u)/] 



where jg is the local Maxwellian distribution 

ft(r,v,t) = n(r,t) 



m 


d/2 


_27rT(r, t) _ 


exp 



m (v — u(r, t)Y 
2T(r,t) 



(22) 



(23) 



This distribution function is parameterized by the nonequilibrium density n , granular temperature T, and flow 
velocity u, which are defined in terms of moments of /: 



dvf, T 



m 
dn 



dvV 2 f, u 



dvvf, 



(24) 
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where V = v — u is the velocity relative to the local flow. The right side of ( |22| ) is a model for the nonlinear Boltzmann 
collision operator. The first term describes collisional relaxation towards the local Maxwellian, with collision rate v. 
The second term describes the dominant collisional cooling effects, where C is the cooling rate. The necessity for this 
term to accurately represent the spectrum of the Boltzmann collision operator is discussed in Ref. ^[ However, it can 
be viewed more simply as an effective "drag" force that produces the same energy loss rate as that produced by the 
inelastic collisions. The parameters v and ( are chosen for good quantitative agreement of the viscosity and cooling 
rate with those obtained from the Boltzmann equation, 

2 (d-i)/2 /T\ 1/2 

imrn* n U) {l - a% ^°' c ' (25) 



1/2 



v Q = k na d - Y [-) , (26) 

where k = 15/16 for d = 3 (spheres) and k = 2 for d = 2 (disks). Both v and £ depend on the density and temperature, 
whose space and time dependence has been left implicit. Further-details motivating these choices can be found in 
Ref. ||, where a somewhat more sophisticated model is described. E3 Equations ( ^2| ) ( ^5|) define the kinetic equation 
to be applied in this work. 

By taking moments with respect to 1, v, and v 2 , this model kinetic equation yields the same form of the macroscopic 
balance equations for mass, energy, and momentum, Eqs. (Q)-(^), as those given from the Boltzmann equation, thus 
confirming the interpretation of £ as the cooling rate. The pressure tensor P and the heat flux q are given by 

dvVVf , q=^ JdvV 2 Vf. (27) 

We now consider the specific problem of steady Couette flow described in the previous Section. The main objective 
is to verify the constitutive equations (p~0|) — (p~3|) from the fundamental definitions (p7|), and to verify the assumed 
constancy for the pressure and shear rate. For the chosen geometry the kinetic equation becomes 

(1 - C* + v y u^ l d y ) f - l -Cd v ■ (V/) = (1 - C)h, (28) 

where 

is a constant that henceforth will characterize the inelasticity dependence. Equation ( |2^ ) allows one to understand 
the structural simplification afforded by the kinetic model relative to the Boltzmann equation. Assuming the hy- 
drodynamic fields of Section O, the local equilibrium distribution function ft is completely specified and Eq. fl28| ) 
becomes a linear inhomogeneous partial first order differential equation that is readily solved with specified boundary 
conditions. This determines explicitly the entire velocity dependence. However, the solution is only formal since the 
assumed fields must satisfy (|24). These are consistency conditions that are needed to justify the assumed forms for 
the hydrodynamic fields. Alternatively, since these fields have been shown to follow from the constitutive equations 
(p^|)-(^3|) for the heat and momentum fluxes in Section [El], it is sufficient to show that these equations are verified by 
use of this formal solution in (|27j) . Both the consistency conditions and the fluxes are verified explicitly in Appendix 



The general solution to Eq. (Pq) is 

/(v,v) = / fl (y,v) + (l-0 / dte-^-i^ e- m h{y,w), (30) 



where T> is the sum of a generator for translations of y and the generator for scale transformation of V, 

V EE Vy^Oy ~ ■ d V . (31) 

The explicit effects of this generator are described in Appendix [A]. The first term of (|30"|), /s(y,v), is a solution 
to the homogeneous kinetic equation obtained from by setting fi — » 0. The detailed form of this contribution 
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is determined by the chosen boundary conditions. The physical boundary conditions are specification of the half 
distributions for velocities directed away from the walls at y — ±Z, given explicitly or implicitly in terms of the 
distributions for velocities directed at the walls. The specific relationship characterizes the motion and temperature 
of the walls. One possibility for Couette flow is diffuse conditions, where the distributions for velocities away from 
the walls are given by a Maxwellian whose parameters specify the temperatures T± and velocities U± of the walls at 
y = ±L. 

In general /s(y,v) has a detailed explicit dependence on the geometry of the system and is responsible for the 
"boundary layer" near the wall. In contrast, the second term of (|30|) is an example of a "normal" solution where all of 
its dependence on y occurs only through the hydrodynamic fields. Here we consider idealized boundary conditions for 
which /s(y, v) — > so that the entire solution is given by the second term of (p0|) . This idealized boundary condition 
eliminates boundary layers and admits the possibility of simple hydrodynamic profiles that are exact throughout the 
system. Furthermore, the solution to the kinetic equation is now normal so that the fluxes become functionals of the 
hydrodynamic fields and (|l|)-(^|) become a closed set of hydrodynamic equations. The idealized boundary conditions 
correspond to the limit T± — > (as illustrated in Fig. [j]) . This is verified in Appendix |b] and will not be discussed 
further here. 

Equation ( |30| ) with these idealized boundary conditions admits an exact solution for steady Couette flow [see Eq. 
(CI) for its explicit form] characterized by the solutions to Eqs. (||), (||), and ( [l5|) given in Section |l|, which are 
rewritten here for the sake of completeness: 



p(y) = n{y)T(y) = const, v Q 1 d y u x = a, (l> 1 d y ) 2 T = -2mPr"/. 



(32) 



The Prandtl number, defined following jl^), is equal to 1 in the kinetic model while Pr = {d— l)/d for the Boltzmann 
equation. The fact that Pr = 1 in the kinetic model is a consequence of the introduction of a single collision 
frequency, that is unable to reproduce simultaneously the exact shear viscosity and thermal conductivity Navier- 
Stokes coefficients. The components of the flow velocity in the y and z directions vanish. The parameter a is the 
dimensionless shear rate and characterizes the relative velocities of the wall. Together with the coefficient of restitution 
a it is a given control parameter in terms of which all transport properties are expressed. The dimensionless parameter 
7 (a, a) characterizes the curvature of the temperature field and is a consequence of both the viscous heating due to 
the shear rate and the collisional dissipation characterized by a. The consistency conditions and forms for the fluxes 
are verified in Appendix [c], confirming that the distribution function ( pfl| ) is an exact solution to the kinetic equation. 
This verification also fixes the functional dependence of 7(0, a) on the control parameters. The details of the analysis 
also are given in Appendix yielding the following implicit equation 



4(1 - C* 
dir 1 / 2 



du t 



dtt 



d-l 



(l + a 2 t 2 )u 2 



(f + 2 7 w 2 ) 1 



(33) 



where 



w(u, t) 



2 i/ 



f 



(34) 



Interestingly, the representation (|33| ) exists only for 7 > or, equivalently, for a equal to or larger than some critical 
value a c of the shear rate. This critical shear rate corresponds to 7 = 0, or uniform temperature (see Eq. ([l9|) 
Setting 7 = in Eq. d33|), one gets 



d C 



2 1 - C* 



(35) 



For the limiting case,a = a c the viscous heating is exactly balanced by collisional cooling and the gas is in a state of 
uniform shear flow £3 This is further discussed at the end of this Section. In the regime of low shear rates and low 
dissipation, Eq. (|3|) yields 7 « (a 2 - d(*/2)/(d + 2) and consequently a 2 « dQ* /2. 

The explicit expressions for the momentum and heat fluxes are also derived in Appendix From them it is possible 
to identify the generalized transport coefficients defined in Eqs. (10)-(|l3|). The viscosity coefficient is 



F v (a,a) = 



4(1 -C* 

ttV2 



du < 



dte~H (l + 2 7 w 2 )" 



(36) 



The viscometric functions characterizing normal stresses are 

2(1 -C) 



*i(a,a) 



du e 



die'* [l - 2 (l 



a t 



(l + 2 7 w 2 



(37) 
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* 2 (o,a) = 



2(1 -C) 
r7 2 7rV 2 



due"" 2 (1 - 2m 2 ) 



(1 + 2 7 w 2 



'0 JO 

The thermal conductivity coefficient is given in terms of 7 and F v by Eq. ( |l6| ) 



Fa 



'/o 



2mAoPr7 



a 2 F„ 



2p 



(38) 



(39) 



Finally, the expression for the cross coefficient <J> in the heat flux is 



$ (a, a) 



8(1-C) 
(d + 2)7rV2 JQ 

2 7W 2 )^ 2 



dtt e 



fit) 



(1 



a 2 * 2 ) 



x (1 + 



(40) 



In summary, the distribution function, hydrodynamic fields, and transport coefficients for the heat and momentum 
fluxes have been determined exactly in terms of the imposed shear rate a and the restitution coefficient a. For small 
a and a the results agree with predictions of the Navier-Stokes hydrodynamics. More generally, they extend the 
description of Couette flow to large spatial inhomogcncity and strong dissipation. The final results are still only 
implicit, but the entire problem has been brought to quadratures. A numerical evaluation of these expressions is 
provided in the next Section for comparison with numerical simulation of the Boltzmann kinetic equation. 

It is easy to check that all results presented in this Section-reduce to those previously derived in the elastic limit 
a = 1 by using the Bhatnagar-Gross-Krook (BGK) modelJatil In addition, ther*€p«lts of this Section include as a 
limiting case those corresponding to the uniform shear flow for an inelastic gas.113113 This happens when the shear 
rate and the inelasticity combine to yield a zero curvature for the temperature profile (7 — ► 0) . Equation ( |35| ) can be 
seen as the relationship between the reduced shear rate and the coefficient of restitution, which are not independent 
in the uniform shear flow problem. The transport coefficients in this limiting case are given in Appendix It is also 
interesting to note that this uniform shear flow state is mathematically equivalent to the corresponding case of elastic 
collisions with an external thermostat force — ^m("V adjusted to control the viscous heating.E£l In this context, Eq. 
( |35| ) (with the adequate change of units) gives the shear rate dependence of the thermostat and this coincides with 
the previous result derived in the uniform shear flow problem for an elastic gas, as expected. 



IV. COMPARISON WITH MONTE CARLO SIMULATIONS 



In this Section we compare the predictions of the kinetic model with computer simulations of the Boltzmann 
equation by means of the Direct Simulation Monte Carlo (DSMC) method. Although this method was originally 
devised for elastic particles]^ its extension to the inelastic case is straightforward.BEj This method has proven to be 
an efficient and reliable tool for solving numerically the Boltzmann equation. 

We have used the adaptation of the DSMC method tehthe steady Couette flow for a granular gas following the same 
steps as those recently worked out in the elastic case,E3 so that the technical details of the simulations are omitted 
here. We have considered a system of hard spheres (d = 3) with three different values for the coefficient of restitution: 
a = 1 (elastic case), a — 0.9, and a — 0.8. For each case, typically four or five different shearing states have been 
taken. Once a steady state is reached, the profiles of the hydrodynamic quantities and of the fluxes are measured in 
the bulk region. From these profiles the transport coefficients are identified. 

The two parameters defining the kinetic model are the collision frequency v and the cooling rate £, given by (pfj). 
These choices were made to provide a quantitative representation of the shear viscosity at the Navier-Stokes level, 
while retaining the relative simplicity of the kinetic equation. To verify the extent to which these choices are valid 
for larger spatial gradients we compare the cooling rate of the kinetic model with that obtained from the Boltzmann 
equation. Both are proportional to (1 — a 2 ). In Fig. ^ we plot the ratio C*/(l — a2 ) obtained from the simulation 
versus the shear rate for different values of a. It is apparent that C* is well approximated by the kinetic model form, 
although the latter tends to underestimate the correct value, especially as the shear rate increases. This supports the 
expectation that the kinetic model retains quantitative as well as qualitative validity even for the extreme conditions 
admitted by the exact solution of the previous Sections. 
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FIG. 2. Plot of the reduced cooling rate divided by 1 — a 2 versus the reduced shear rate a, as obtained from simulation 
for a = 0.98 (triangle), 0.95 (diamond), 0.9 (circles), and 0.8 (squares). The dotted line represents the local equilibrium 
approximation C*/(l — °^ 2 ) ~ 5/12 used in the kinetic model. 

The simulation results indicate that the pressure p is practically constant, in agreement with the first assumption in 
([32|). Also, the actual shear rate d y u x multiplied by T 1 / 2 is practically constant, allowing identification of a reduced 
shear rate a in agreement with the second equality in (|32|). We have also checked that the temperature profile is well 
characterized by the third equation in ( |j2|) (with Pr = 2/3). The parameter 7 characterizing the curvature of the 
temperature profile is plotted in Fig. || as a function of the shear rate for a — 0.8, 0.9, and 1. The agreement between 
theory and simulation is quite good. Surprisingly, this is true even for strong inelasticity. There is a remarkable 
influence of the degree of dissipation on 7 at a given value of the shear rate. As the coefficient of restitution a 
decreases, so does the temperature variation across the system, so tjbat-jthc state tends to resemble the uniform shear 
flow, in which case the kinetic model is known to be quite accurate. Eata Figure]^ also confirms that, for a given value 
of a, there exists a critical shear rate a c such that 7 — > + as a tends to a c from above. For shear rates smaller than 
the critical value a c the collisional cooling effect dominates over the viscous heating, and the temperature presents a 
minimum instead of a maximum, which could be characterized by a negative 7.r 9 H 30 l However, the solution obtained 
here does not extend to negative 7, as noted above. The dependence of a c on the coefficient of restitution a is shown 
in Fig. ^. The simple prediction ( p5| ) reproduces quite well the simulation data. Figure || can also be interpreted as 
the two dimensional parameter space of the inelastic Couette flow. The elastic Couette problem corresponds to the 
line a = 1 and arbitrary a. On the other hand, the uniform shear flow for a dilute granular gas is represented by 
the line a = a c {a) for arbitrary a. The case studied in this paper applies to the region enclosed between both lines, 
namely, for arbitrary a and a, provided that a > a c (a). 
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Now we compare the simulation results and the theoretical predictions for the five generalized transport coefficients 
in the domain a > a c . For the sake of clarity, we have also plotted the curves representing the transport coefficients 
for the uniform shear flow when eliminating the coefficient of restitution in favor of the shear rate. In Fig. we 
plot the shear rate dependence of the viscosity function F v defined by Eq. (10). Needless to say this is the most 
relevant quantity in a shearing state. Regardless of the value of a, shear thinning effects are present, i.e., F„ decreases 
as the shear rate increases. This rheological behavior is quantitatively well described by the kinetic model, again 
especially for large dissipation. In general, however, the model tends to overestimate the value of F v . The two 
viscometric functions (jll) are plotted in Figs. ^|and|j. While the agreement between the theoretical predictions and 
the simulation results in the case of ^1 is similar to the one observed with F v , the agreement in the case of VP 2 is only 
qualitative. In particular, the model succeeds in capturing the non-monotonic behavior of \&2 in the inelastic case. 
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FIG. 5. Plot of the viscosity function F v as a function of a 2 , as obtained from the kinetic model (lines) and from simulation 
(symbols), for a — 1 (solid line and triangles), 0.9 (dashed line and circles), and 0.8 (dotted line and squares). The dash-dotted 
line refers to the corresponding transport coefficient on the critical line (uniform shear flow). 
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FIG. 7. The same as in Fig. pmut for the second viscometric function \p2- Note that $2 = on the critical line. 

The transport of energy is characterized by the coefficients F\ and $, defined by Eq. (|l3|). They are plotted in 
Figs. |^ and ^, respectively. As happens with F n , the generalized thermal conductivity F\ decreases as the shear 
rate increases. The accuracy of the kinetic model predictions for this quantity is less satisfactory than in the case 
of the viscosity function. This is a consequence of the fact that the model slightly underestimates the cooling rate 
(Fig ||) and the curvature parameter 7 (Fig. ||), while it slightly overestimates F v (Fig. ||) and so, according to 
Eq. (|39|) , all these discrepancies contribute to magnify the inaccuracy of the value of Fx predicted by the model as 
well. In the case of the coefficient $, the agreement between theory and simulation is quite good. This is rather 
satisfactory if one takes into account that this is a coefficient measuring complex coupling effects between the velocity 
and temperature gradients, which are absent at the Navier-Stokes regime. Also note that the kinetic model yields 
$(a = 0, a = 1) = —2.8, while $(a = 0, a = 1) = —3.5 in the Boltzmann equation. This is again due to the fact that 
relaxation effects are described by a single frequency in the BGK model. 
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FIG. 8. The same as in Fig. B but for the thermal conductivity function F\. 
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FIG. 9. The same as in Fig. H but for the cross coefficient <E>. 



V. DISCUSSION 



The exact macroscopic balance equations for mass, energy, and momentum become closed hydrodynamic equations 
when the energy and momentum fluxes are specified in terms of the hydrodynamic fields. In Section || specific forms 
for these fluxes were postulated to hold for Couette flow, and the resulting hydrodynamic equations were solved 
exactly for the temperature, pressure, and flow velocity. In Section [II an exact solution to a model kinetic equation 
was obtained for Couette flow and the assumed forms for the fluxes was verified. In addition, it is shown in Appendix 
^ that the hydrodynamic fields obtained directly from this distribution function are the same as those obtained 
from solution to the hydrodynamic equations. Thus, both hydrodynamic and kinetic descriptions have been verified 
as exact and self-consistent. The kinetic theory provides the distribution function in phase space, Eq. (CI), as an 
explicit functional of the hydrodynamic fields and therefore allows in addition the calculation of non-hydrodynamic 
properties in terms of them. These fields are given by Eqs. ( |l9| ) and with the parameter 7 (a, a) determined by 
the consistency condition (|33|). 

Some comments to summarize these results and provide perspective are as follows: 

• The model kinetic equation considered is a generalization of the familiar BGK equation used to model the 
Boltzmann equation for elastic collisions, and reduces to it for a = 1. The utility of the BGK equation to 
address complex states not accessible via the Boltzmann equation is well-established. The results here and 
those of Refs. |lj and [H] confirm its value for granular media as well. Although these kinetic model equations 
are structurally simple, they are highly nonlinear due to the implicit dependence of fe on the distribution 
function via the hydrodynamic fields. 

• The kinetic model equation provides a coupled set of singular nonlinear integral equations for the hydrodynamic 
fields. Incorporation of specific boundary conditions is straightforward and an advantage of this description. 
However, explicit results generally require approximations near the homogeneous state or numerical solution. 
Exact analytic results valid for large spatial gradients and arbitrary inelasticity are rare even for the kinetic model 
equations. They require idealized boundary conditions such that boundary layer complications are avoided and 
"simple" profiles are possible in the appropriate variables. 

• The idealized boundary conditions considered here are homogeneous, such that the solution is "normal" — all 
space dependence occurs entirely through the hydrodynamic fields. Since the fluxes in macroscopic balance 
equations (0)-(H) are determined from moments of this distribution, they can be calculated as functions of the 
hydrodynamic fields. Use of these in the macroscopic balance equations then gives a closed set of equations for 
the fields. The analysis here provides a non trivial example of the existence of hydrodynamics for a strongly inho- 
mogeneous state together with the solution to those equations. The solution is valid even for strong dissipation 
subject to the condition 7 > or a > a c = \/ d(* (a) /2[1 — C*(°0]- 
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• The results, derived here include as a limiting case those previously obtained for the uniform shear flow 
problem£jll3 This corresponds to an exact balance between the viscous heating and the inelastic cooling, 
resulting in a uniform temperature (7 — > + ). From that point of view, the critical shear rate a c gives the 
relationship between the shear rate and the coefficient of restitution in the uniform shear flow. 

• The existence of a critical shear rate a c below which the analysis does not apply restricts the conditions such 
that viscous heating dominates collisional cooling. The results do not apply, therefore, to the simple-ease of zero 
shear rate and constant temperature walls. In this latter case, simulations and asymptotic analysisoEJ suggest 
that the temperature profile is very close to that of (|3^) but with 7 < 0. The proof for this is still not available. 

• Although the kinetic model is only a crude representation of the Boltzmann equation, it does preserve the most 
important features for transport, such as the reference homogeneous state and the macroscopic conservation 
laws. Otherwise there are only two adjustable parameters, the collision rate v and the cooling rate (, to improve 
correspondence with the Boltzmann equation. Here they have been chosen to fit the Navier-Stokes order viscosity. 
Surprisingly, this is sufficient to allow good agreement, with the Boltzmann results for rheological properties far 
from equilibrium. This is true for elastic collisionsOES and now has been confirmed here for inelastic collisions. 
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APPENDIX A: ACTION OF THE EXPONENTIAL OPERATOR 

In this Appendix we evaluate the action of the operator e~ vt , where V is given by Eq. ([H]). To do this, it is 
convenient to refer the velocities of the particles to a local Lagrangian frame V = v u and to introduce a scaled 
space variable as 

V dy'My')- (Al) 

yo 

In terms of this variable, the velocity and temperature profiles become simply 

u x (s) = as, T (s) = T - 7717s 2 , (A2) 

where we have taken into account that Ao = id + 2)r)o/2m in the BGK kinetic model. In the new Lagrangian frame 
the operator T> becomes 

V = V y d s - aV v dv x - yV • dv, (A3) 

where the derivative d s is taken at constant V. 

Let us consider the operators A = —V y d s) B = \C,*~V • &v, and C = aV y dv x - It is easy to see that the operators 
A + B and C commute, namely, 

[A + B,C] = (A + B)C- C(A + B) = 0. (A4) 
Thus, the action of the operator e~ vt = et^+^+c)* on a given function F(s,\) is 

e (A+B+C)t F ^ y) = e (A+B)t e Ct F = e Ct e (A+B)t F 

= e (A+B)t F(s,Y + tz ■ V), (A5) 

where a is the matrix with elements Q>ij — CL$ix$jy ■ It remains to evaluate the action of the operator + " on an 
arbitrary function of the velocity, h(s, V): 
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e( A+B ^h(s,V) = e m H(t,s,V). 
An equation for H(t, s, V) follows by differentiating both sides to get 

d t H(t, s, V) - e- Bt Ae Bt H{t, s, V) = 0. 
Using the property for scale transformation 



(A7) simplifies to 

d t H(t,s,V) + e-*r*V y d a H(t,8,V) = 0. 
Integrating this with the initial condition H(0, s, V) = h(s, V) gives 



H(t,s,V) = exp 



|-(l-e-^ C *')^9 s 



/i(s,V). 



Use of this in (A6) gives the desired result 

e {A+B)t h{s, V) = exp ( l(*tV • <9v ) exp 



v 2 

exp [— r 



1 - e-*<** ) V y d s 



h(s,V) 



where 



{t)V y ds\ h{s, e5 C** V ) = ft ( s - T(t)V B , e' c "*v) 



(t) = 1 ( e*<*« - 1 



Finally combining Eqs. ( Rq ) and (All) gives the result 

e - ra F(s, V) = e (A+B+c)t .F(«, V) = F ( s - r{t)V y , e* c ** (V + ta ■ V) 



(A6) 
(A7) 

(A8) 
(A9) 
(A10) 



(All) 



(A12) 



(A13) 



APPENDIX B: IDEALIZED BOUNDARY CONDITIONS 

The boundary conditions for the general solution (KG) are imposed by the choice of fs, which is a solution to 



or, equivalently, 



i-C --C +v) f B = o 



C* 



dsfs = - V 1 - C* - tjC* - aV w ^ - ^-V • dv) f B , 



where we have made use of Eq. ( A3 ) . This can be integrated to get 
/b(s,V) = e(V r „)exp 



S + "' 'l-C-iC-aVydy.-^V.dv 



+6(-K)exp 



$_(V) 
*+(V). 



(Bl) 



(B2) 



(B3) 



Here 0(x) is the Heaviside step function, s — ±£± correspond to the location of the walls at y = ±L, and $±(V) 
are arbitrary functions. These functions are fixed by the boundary conditions at s — ±£±: 



Q(±V y )f(TC T ,V) = Q(±V y )$ T (V). 



(B4) 
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Idealized boundary conditions are constructed as follows. First we write the boundary conditions as 

e(±V y )f( T C T ,V) = e(±V v )^(m/2T T ) d /^ T (V/V27VM) , (B5) 

where 

n T ee (2T T /my 1/2 J dV e( T V y )\V y \f(T^, V), (B6) 

C t ee|^O(±^)|4|^(0, (B7) 
V=f(£) being the distributions of velocities away from the walls normalized as 

The usual choice is the Maxwell-Boltzmann distribution = <£>-(£) = r K~ d l 2 cxp(— £ 2 ). Now consider the limit 

of infinitely cold walls, T± — > 0. In that limit $±(V) — > since <p±(£) must vanish if |£| — » oo. Consequently, 

lim / B (s,V) = 0. (B9) 

APPENDIX C: CONSISTENCY CONDITIONS AND CALCULATION OF THE FLUXES 

1. Generating function 

The consistency conditions for the assumed forms of the hydrodynamic fields are the d + 2 equations of (|24|) . 
Evaluation of the right hand sides of these equations requires the explicit form for the distribution function in terms 
of these fields. Using the results of Appendix [a|, the solution ( |30| ) with the idealized boundary conditions discussed 
above becomes 



/(s, V) = (1-C) / dt e -(W*-!0* e atV,8v„ e -r(t)V,fl. / £ ( s , e iC** V ) 
Jo 

/■oo 

= (1 - O J dt e-^'-i ^f e (a - V y r(t), e^* 4 (V + ta • V)) , (CI) 
where the local equilibrium distribution is 

Ma,V) = ^e(T( s ))v- {d+2 \ s )7r^ 2 exp[- (V/vo(s)) 2 } (C2) 



in 



with «§(s) = 2T(s)/m. 

The velocity integrals for the consistency conditions and for the fluxes are all low order moments of this distribution 
function. They can be obtained from appropriate derivatives of the generating function 

G(s,k) = J dVe ik V /( S ,V). (C3) 

Making use of the solution (CI), the generating function becomes 

G(s,k) = (1-C) J™ dte-^S'-i W J dVe lk - v e atv »^e^W y ^ f e (s,ei <'*V) 

= (1 - C) f°° dt e^ 1 -^* f dVe 4k ' (t) - v e- Tl(t)v ^ fe(s,V), (C4) 



o 



where 
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k'(t) = e -H** e - atfe ^^k, n (<) = (i - e - * 

In the last step of Eq. ( |C4| ) we have made the change V — > *V and have used the general property 

J dV F 1 (V)e atv y dv - F 2 (V) = J dV F 2 (V)e~ atv y dv - F^V). 
Integration over V_l = V V y y yields 

G(s,k) = G+(s,k) + G_(s,k), 

where 

G±(s,k) = ^(1 - C) / dV„ / iife-( 1 -f , ) t e ±< *iW l, «e^W v ' 9 'e- , ?^WF(« 1 k / i (t)) 1 
m Jo Jo 

with 

F( S ,kl(0) = e(« 2 ( S ))V 3 (s) 7 r- 1 / 2 exp 



-J*± a (*)«oW 



(C5) 



(C6) 



(C7) 



(C8) 



(C9) 



We will focus on G+(s,k) since G_(s,k) can be evaluated as G_(s,k) = G + (— s, — k). 
Next, define the change of variables in the t integration 

z = V v Ti{t), dz = [Vy - -4r- I dt, t = -i- In [ 1 - C 



2K 



(CIO) 



The function G+(s, k) becomes 



2p 



G + ( s ,k) = ^(i-o / 

xe- v S^ s -^F{s ~ z,k' ± (t)) 

= due-^u- 1 dze ~(^C)tjk' y (t)u V0 (s-z) 

>o Jo 

c* 



m 

x6 I 1 - 



1 - 



2 uvq{s — z) 



F(s-z,k' ± (t)). 



2 uvq(s — z) 

A change of variables in the V y integral has been made, V y — > uvq(s — z). Accordingly, the variable t becomes 



t= In 

c* 



1 - 



C*z 



2uvq(s — z) 



Next, for the z integral change variables to 



v (s - z) ' 



(Cll) 



(C12) 



(C13) 



As a consequence, z as a function of w is z = z+(s,w), where z+(s,w) is the positive root of the quadratic equation 
obtained by squaring QC13 ): 



It can be easily shown that 



The generating function is now 



(l + 2~/w 2 ) z 2 - 4'jw 2 sz - w 2 v 2 (s) = 



£ = ^)(^) 3 [i-^ s inr( s ) 



(C14) 



(C15) 
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G+(»,k) 



j(a)7rV2 



70 

l-^0,lnT(a)' 



2 it 



2 u 



-fei 2 (t)( j:+ /2«;) 2 



where t is now a function of u and 



-Alnfl-^ 



Finally, change variables in the w integration to t, so that 



2u 



This leads to the explicit result 



G+(s,k) 



die~ (1 ~ r) *A+(s,k,u,i), 



where 



A+(s, k, u, i) 



z+(s,u;) 



d s lnT(s) 



exp 



z+(s. w) k\ Z+(s, w) 
~T(t) 4" u 2 r 2 (t) 



It must be recalled that w is a function of u and t given by (C18) 
The function G_(s,k) = G+(— s, — k) is then 



G_(»,k) 



oo 

-"" ' dte-( 1 - c *)*A_(s,k,u,t), 
o 



with 



A- (s, k, u, t) 



z_(s, w) 



cUnT(s) 



exp 



Z-(s, w) k\ z_ (s, w) 



T(t) 4 U 2 T 2 (t) 



where Z-(s,w) = —z + (~s,w) is the negative root of the quadratic equation (C14). Real roots require 

(4 7 i 

which is satisfied for positive T(s) and positive 7. 



sw 2 ) 2 + 4 (1 + 2jw 2 ) 2w 2 ~^- > 0, 



2. Consistency conditions 

The consistency condition for the density becomes 



n(s) = /" dvf = G(s, k = 0) 



l-^9 s lnT( S ) 



T( S )- 

Use has been made of the identity 



l-—d s \nT(s) 



= 2, 



l-^-5 s lnT(s) 
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which follows from the explicit forms for the roots of ( |C14 ). The result (C24) is consistent with the required equation 
of state defining p. 

The consistency condition for the temperature is 



dn(s)T(s) =mj dvV 2 f = -m [%G(s, k)] k=0 . 



Direct evaluation gives 
1 = 



m 
dT(s 

x < zl 



due~ u ' I dte- f 



d-l 



+ {l + a 2 t 2 ) 



+ z 2 _ 



l-—d s \nT(s) 



(C26) 



(C27) 



Using the identity 



l- Z -±d s \nT{s) 



l-— ds \nT( S ) 



4T(s) w 2 
m 1 + 27W 2 



one gets the hnal result 



1 



dn 



due u die 
[i Jo 



d-l 



+ (l + a 2 t 2 ) u 2 



(l + 27W 



2\"1 



This condition is enforced by using it as the definition of 7. 

Finally, the consistency conditions for the flow velocity components arc 



= mf dvV i f=[d kt G(s,k)] 1e=0 . 



(C28) 



(C29) 



(C30) 



This is clearly satisfied for the z component. The x and y components are proportional to integrals including the 
term 



l-^cUnT(s) 



+ -— 



1-— d s \nT(s) 



(C31) 



which identically vanishes. 

This co mp letes confirmation of the consistency conditions for the hydrodynamic fields. In summary, the distribution 
function ( |Cl| ) with the hydrodynamic fields ( |32|) constitutes an exact solution to the kinetic equation. 



3. Calculation of the fluxes 



The macroscopic transport properties of the steady Couette flow are given by the momentum and heat fluxes, Eq. 
(p7j) . In terms of the generating function these are 



[9 fci a fcj G( S ,k)] k=0 , q(s) 



iim[^G( S) k)] k=0 . 



(C32) 



Following similar steps as in the consistency condition for the temperature, the nonzero elements of the pressure 
tensor can be evaluated. They are given by 



Pr 



2p 



TTa^-O / due 



dti 



1 + 2 (atu) 2 } (1 + 2 7 w 2 ) 1 , 



(C33) 



dt e~* (1 + 2jw 2 y 



(C34) 



2p 



du e 



dt e~ f (l + 2~/w 2 ) 1 



(C35) 



19 



dte'H (l + 2jw 2 



(C36) 



It is readily verified that these results satisfy the relationship P — \ [P X x + Pyy + {d — 2)P ZZ \. 

Let us consider now the heat flux vector. Note that, although the temperature gradient is only directed along 
the y direction, the presence of the shear flow induces a nonzero x component of the heat flux. Thus the nonzero 
components are q x and q y . Using ( C32| ), q y is found to be 



q y {s) 



mp 



2T( S )7T 1 /2 



POO POO 

— (1-C) / due~ u \ dte-^+^'^w- 3 



d-1 



+ {l + a 2 t 2 ) 



l- Z -±d s \nT(s) 



+ z°_ 



87PS 

where use has been made of the identity 



1-— d s hxT(s) 



d-1 



+ {l + a 2 t 2 ) 



2^ u 2 



w(i + 2-fw 2 y 



1 - -fd s lnT(s) 



l-—d s \nT(s) 



I67S 



T(s) 



w 



m (1 + 27W 2 ) 2 



This can be written as 



9v = -^i72( 1 -^ ) J Q due 
dp(* 4a 2 ps 



dte~ l 



d-1 



+ (l + a 2 t 2 )u 2 



Of 



1 + 2 -yw 2 



+ ^7r(l-C*)/ due- u 2 I dte-H(l + 2 7 w 2 ) 



(C37) 



(C38) 



(C39) 



An integration by parts in the t integral and the consistency condition (C29) have been used to obtain the last equality. 
Comparison of this result with that for P xy above shows 



1y = 7, — 7, 1" aP xy O s T. 



2m7 V 2 



(C40) 



This result is in fact required by the hydrodynamic equations to support the forms for the hydrodynamic fields. 
In a similar way q x {s) is calculated, 



9x(s) 



amp 



2T( S )7T 1 /2 



p 00 poo 

— (1-C*)/ due- u2 u dte-^+^'hw' 3 



d+1 



+ (l + a 2 t 2 ) u 2 



x < z 



l- 7 ^d s \nT{s) 



+ z-_ 



l-—d s \nT( S ) 



4p 



TflTT 



™a(l-C*) / due- u2 u dt e-( 1+ ^')hw 



d+1 



+ (1 + a 2 t 2 ) u 2 



x (l + 2 7W 2 ) 2 <9 S T. 



(C41) 



APPENDIX D: TRANSPORT COEFFICIENTS IN THE UNIFORM SHEAR FLOW LIMIT 



In this Appendix we derive the explicit expressions for the transport coefficients along the critical shear rate a c (a). 
They are obtained by taking the limit 7 — > + in the corresponding expressions of Sec. III. Equations (36)-(|40|) 
simply yield 



F v (a c ,a) = 1-C = 



d + 2aV 



(Dl) 



*i(o c ,a) = -2(l-C*) = 



2d 



d + 2aV 



(D2) 
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* 2 (a c ,a) = 0, 



(D3) 



$>(a c ,a) 



2 (1-C)(4 + 3Q 
d+2 (1 + C*) 4 (2 + C*) 4 



7(2 + 3C* + C* 2 ) 2 + I8a 2 c (8 + 12C* + 5C 



*2x 



(D4) 



In order to obtain F\ from Eq. (f39|), we need to evaluate a 2 and F v to first order in 7. From Eqs. (|3^) and ( |36| ) one 

gets 



2 1-C 



A(Ch 



(D5) 



F, = l-C*-B(C*)7 + -"» 



(D6) 



where 



A(C) = 



(1 + C) 3 (2 + C*) 3 



2x2 



(d + 2)(2 + 3C*+C* ) 



-6^(24 + 48C* + 33C* 2 + 9C* 3 + C* 4 ) 



(D7) 



g(n = 12 d + m2+c) 2 (6+6C+r 



Thus, the thermal conductivity is 



fj(o c ,a) 



(l-C)A(C)-agg(C) 
d+2 



(D8) 



(D9) 



Equations (D1)-(D3) coincide with those previously derived in Ref. [12| However, E qs. (D4) and (|D9j) are new 
results. Despite the fact that there is no heat flux in the uniform shear flow, Eqs. (D4) and (DE) are intrinsic 
transport coefficients characterizing the state of the system. 



1 C. S. Campbell, Ann. Rev. Fluid Mech.:22, 57 (1990). 

2 J. J. Brey, M. J. Ruiz-Montero, and D. Cubero, Phys. Rev. E:54, 3664 (1996); J. J. Brey, D. Cubero, and M. J. Ruiz-Montero, 
ibid. :59, 1256 (1999); S. Luding, M. Miiller, and S. McNamara, in World Congress on Pa rticle Technology (B righton, 1998, 



CD: ISBN 0-85295-401-9); J. M. Montanero and A. Santos, Gran. Matt.:2, 53 (2000) and |sond- mat/0002323 . 

3 J. J. Brey, J. W. Dufty, C. S. Kim, and A. Santos, Phys. Rev. E:58, 4638 (1998). 

4 N. Sela and I. Goldhirsch, J. Fluid Mech.:361, 41 (1998). 

5 V. Garzo and J. W. Dufty, Phys. Rev. £:59, 5895 (1999). 

6 J. J. Brey, J. W. Dufty, and A. Santos, J. Stat. Phys.:97, 281 (1999). 
7 R. Zwanzig, J. Chem. Phys.:71, 4416 (1979). 

8 A. Santos, J J. Brey, and V. Garzo, Phys. Rev. A:34, 5047 (1986). 

9 J. J. Brey, A. Santos, and J. W. Dufty, Phys. Rev. A.36, 2842 (1987). 

10 C. S. Kim, J. W. Dufty, A. Santos, and J. J. Brey, Phys. Rev. AAO, 7165 (1989). 

11 J. M. Montanero and V. Garzo, Phys. Rev. E:58, 1836 (1998). 

12 J. J. Brey, M. J. Ruiz-Montero, and F. Moreno, Phys. Rev. #55, 2846 (1997). 

13 J. M. Montanero, V. Garzo, A. Santos, and J. J. Brey, J. Fluid Mech.:38Q, 391 (1999). 

14 J. T. Jenkins and S. B. Savage, J. Fluid Mech.:130 , 187 (1983); J. T. Jenkins and M. W. Richman, Phys. Fluids:28, 3485 
(1985); J. Fluid Mech.:192, 313 (1988). 

15 C. K. K. Lun, S. B. Savage, D. J. Jeffrey, and N. Chepurniy, J. Fluid Mech.:140, 223 (1984). 

16 M. A. Hopkins and H. H. Shen, J. Fluid Mech.:244, 477 (1992). 

17 I. Goldhirsch and M. L. Tan, Phys. Fluids:8, 1753 (1996); N. Sela, I. Goldhirsch, and S. H. Noskowicz, ibid.:8, 2337 (1996). 

18 C. Cercignani, "Shear flow of a granular material," J. Stat. Phys., to be published. 



21 



19 M. W. Richman and C. S. Chou, J. Appl. Math. Phys.:39, 885 (1988). 

20 T. N. Hanes, J. T. Jenkins, and M. W. Richman, J. Appl. Mech.:55, 969 (1988). 

21 C. K. K. Lun, Phys. Fluids:8, 2868 (1996). 

22 M. Babic, Phys. Fluids:Q, 2486 (1997). 

23 S. Chapman and T. G. Cowling, The Mathematical Theory of Nonuniform Gases (Cambridge University Press, Cambridge, 
1970). 

As an illustrative example, let us consider a system of bronze spheres of diameter a = 0.35 mm (m = 2 x 10~ r kg) with a 
solid fraction <f> = 3.2%, and a pressure p = 3 x 10 3 Pa, recently used in a low gravity experiment [E. Falcon et al, Phys. 
Rev. Lett. :83, 440 (1999)]. This implies v = 0.94 kHz. 

25 The model kinetic equation of Ref. fi uses the exact local homogeneous cooling state of the Boltzmann equation, for which ft 
is known to be an excellent approximation. In addition, it uses the exact cooling rate of the Boltzmann equation, which is a 
bilinear functional of /. If this functional is evaluated using ft, then the form for ( used here is obtained. Again, this appears 
to be a good approximation. Confirmation of the quantitative agreement between the kinetic model and the Boltzmann 
equation for the state considered here is given in Section [tv| . 

26 A. Santos and V. Garzo, Physica A.21S, 409 (1995); V. Garzo and A. Santos, ibtd.:213, 426 (1995). 

27 G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows (C larendon, Oxford, 1 994). 

28 J. M. Montanero, A. Santos, and V. Garzo, Phys. Fluids:12, 3060 (2000) and |sond-mat/0003364 J. M. Montanero, M. 
Alaoui, A. Santos, and V. Garzo, Phys. Rev. EAQ, 367 (1994). 

29 J. J. Brey and D. Cubero, Phys. Rev. E:57, 2019 (1998). 

30 R. Ramirez, D. Risso, R. Soto, and P. Cordero, Phys. Rev. E:62, 2521 (2000). 



22 



